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In models of coupled dark energy, in which a dark energy scalar field couples to other matter 
components, it is natural to expect a coupling to the inflaton as well. We explore the consequences 
of such a coupling in the context of single field slow-roll inflation. Assuming an exponential potential 
for the quintessence field we show that the coupling to the inflaton causes the quintessence field to 
be attracted towards the minimum of the effective potential. If the coupling is large enough, the 
field is heavy and is located at the minimum. We show how this affects the expansion rate and the 
' slow-roll of the inflaton field, and therefore the primordial perturbations generated during inflation. 

, We further show that the coupling has an important impact on the processes of reheating and 

preheating. 
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I. INTRODUCTION 



The nature and origins of dark energy, the energy component which is responsible for the observed accelerated 
expansion of the universe, remain a mystery. Although the observations can be accounted for by the cosmological 
constant .scalar fields [l|, l2| and modified gravity theories |3fl have also been suggested (for reviews and references 
see e.g. [3, @, One of the major aims of modern cosmology is to determine the properties of dark energy, many 
of the parameters of which (equation of state, matter couplings etc.) can be constrained by considering data from 
^jjy supernovae at high redshifts, observations of anisotropics in the cosmic microwave background radiation (CMB) and 
the large scale structures (LSS) in the universe. At even higher redshifts and in the early universe, constraints from 
varying constants and big bang nucleosynthesis (BBN) give further constraints on the dark energy evolution 0, [1[ - 
Some models of dark energy can even be tested in the laboratory 0, [13, EH, EH EH- I n this paper, we study the 
■ impact of couplings of a quintessence-like scalar field to the inflaton field, the scalar field responsible for an accelerated 
expansion in the very early universe. In this important epoch, the seeds for the structures we observe in the universe 
were created. It is usually assumed that dark energy is not important during inflation. In the case in which dark 
energy and the inflaton field are not coupled, the vacuum expectation value (VEV) of the dark energy field is driven 
£N| ■ by q uantum fluctuations to large field values, but otherwise there are no consequences for the inflationary dynamics 
y—i ' ([3> see a l so Oil)- We show that this is not necessarily the case if the inflaton field couples to the dark energy field. 

In models such as coup led quintessence ([3, E3|; see Q for a recent review) or quintessence models with a growing 
matter component [18| . dark energy couples to at least one species, which is thought to be the decay product of the 
inflaton field. Therefore, it is natural in these types of models to consider a coupling between dark energy and the 
inflaton field as well. We will show in this paper that for large enough couplings (to be specified below) one can expect 
modifications to the predictions of the spectral index, its running and the tensor-to-scalar ratio. We find also that 
the details of the physics of reheating and preheating are affected by the presence of a coupling between dark energy 
and the inflaton. The paper is organized as follows: in Section |TT] we describe our model and study the inflationary 
epoch and discuss the effect of the quintessence field on the primordial perturbations. In Section ITTT1 we discuss the 
consequences for reheating and preheating. We conclude in Section HVl 



II. SLOW-ROLL INFLATION IN THE PRESENCE OF COUPLED DARK ENERGY 

The theory we consider in this paper is specified by the action 

r /vf2 

with 



1 

A 2 (Q) 



£q = -^"(dnQ)(d v Q)-V(Q) 



£0 = - AiQYU^) 



2 



Here, 1Z is the Ricci scalar, <p is the inflaton field, Q is another scalar field, possibly playing the role of dark energy 
and g is the determinant of the metric tensor. The appearance of the coupling function A(Q) follows directly from the 
fact that we focus on a scalar tensor theory where matter, i.e. the inflaton field here, couples to the rescaled metric 
9nv = A 2 (Q)g filJ . While these equations are valid for any potential U(<p), we will use as an example the standard 
chaotic inflation potential 

U{<j>) = \m 2 <f. (1) 

For the coupling function we choose A(Q) — exp(/3Q/Mpi), (as in, for example, [H, E3j E^l ) • After a field redefinition, 
the (effective) mass of the inflaton field is nothing but A(Q)m, which grows as Q rolls down along the potential 
towards large values. Additionally, we choose the potential for the quintessence field to be an exponential potential, 
i.e. 

V(Q) = M P1 exp(-XQ/M Pl ), (2) 

with A positive. 



A. Slow-Roll Inflation 



Let us now consider the inflationary period. Since Q couples to <fr, the inflationary dynamics will be modified by the 
coupling. The coupling can potentially ruin inflation and therefore we first consider the conditions on the parameters 
in order to obtain a period of slow-roll inflation. 

The equations of motion for the fields <fr and Q in a homogeneous and isotropic universe are 

BV 

Q + 3HQ+^=A 2 p( ( f> 2 -4A 2 U), (3) 
dU ■ ■ 

</> + 3H<j> + A 2 — = -2/3Q0. (4) 



h2 = ^ttt (^ 2 + lQ 2 + v + A i u). (5) 



The Friedmann equation reads 

3M 2 X ^ 2 T 2 

The Q'fidcl moves in an effective potential, which is given by the bare potential and a part coming from the coupling 
to if). Taking A positive, the sign of determines whether the effective potential has a minimum or not. For positive 
j3, a local minimum exists, whereas there is not one in the case of a negative (3. If the minimum does not exist, the 
effective potential is of a runaway form and the discussion of the field behaviour will be similar to that in [T3|. In 
this paper, we will assume the existence of a minimum, that is, we consider the case of a positive (3. Assuming then 
slow-roll of the scalar field cp, this determines the value of Q at the minimum of the effective potential to be 

gmin = 473TA ln (wS))' (6) 

together with the condition for the minimum 

XV(Q miu ) = A13A A U(£). (7) 
For later use, we also state A(Q m { n ) during slow-roll inflation: 

/ x Af 4 n/3/(4/3+A) 

A ^> - (w§y) ■ < 8 » 

During inflation, in which we assume that both of the fields are rolling slowly, using eqn. (JT]) the Friedmann equation 
takes the form: 
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We will justify below the assumption that both fields roll slowly. Note that the Q-field contributes to the expansion 
rate with an amount depending on /3 and A. The mass of the Q-field can be found to be 



d 2 V 
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In the case of a slowly-rolling (/>-field, this gives 
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where in the last line we have used the Friedmann equation. For large enough (3, the field is rather heavy (m 2 s > H 2 ) 
and will therefore settle into the minimum of the effective potential. To get a rough idea on how small (3 can be, we 
can solve the equation m 2 s /H 2 = 1 and obtain /3 CI it = (we remind the reader that (3 has to be positive for a 
minimum to exist, since we are assuming A positive). This equation implies that, even when [3 is rather small, the 
Q-ficld sits at the minimum of the effective potential. We find numerically that Q settles into the minimum even for 
(3 as small as 0.05 (see section Hi CI) . 

It is possible to obtain a degree of analytical insight into the behaviour of the system during the inflationary period 
by studying the slow-roll regime and deriving the slow-roll parameters. Firstly, we show that the extra friction term 
in eqn. @ is negligible during inflation. Throughout this analysis, we will assume that the value of /3 is such that Q 
settles into the minimum of the effective potential. Thus one can write Q = Q m in- Using eqn. ([6]) we find (' = d/d(j)) 
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where eqn. JT]) was used in the last step. Therefore we find 
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In general, we see that the extra damping term is proportional to </> 2 and we can therefore neglect it: 
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This means that the condition for slow-roll differs from the standard case by a factor of A 2 : 
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The first slow-roll parameter is defined by e 
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jp- and in terms of the fields can be found to be 
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where eqns. (TTBl and ([9]) were used. We can find the rj parameter in a similar way. Differentiating eqn. (fTSf yields 
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The requirement that |0| <C \H(f>\ gives 
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B. Cosmological Perturbations 



We now derive the relevant equations for the cosmological perturbations. We show first that the perturbations in 
the Q-field are much smaller than the perturbations in the 0-field. This is because the Q-field is a heavy field (i.e. 
m e fi > H) and it is well known that for these fields perturbations are suppressed. The fluctuations of a massive 
scalar field during inflation satisfy an equation composed of oscillatory and non-oscillatory parts [2(| • The expansion 
of the universe causes the wavelength of the perturbations to be streched, so that the former can be neglected as the 
average contribution of oscillations averages to zero. The amplitude of the resulting power spectrum is suppressed 
by (H/m e g) and decreases rapidly for large wavelengths. For completeness, however, we will study numerically the 
perturbations for both fields. The inflaton field is light (m < H) and its perturbations are not suppressed. 

In the longitudinal gauge and in the absense of anisotropic stress, the scalar perturbations of the FRW metric can 
be expressed as (see e.g. |2]| ) 

ds 2 = -(I + 2^)dt 2 + a 2 (I - 2^)S ij dx i dx j . (20) 



The perturbed Einstein equations give 
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where 



3 tf(* + **) + _* = -^p, (21) 

* + 2T* = -^, (22) 

* + 4# * + (3H 2 + 2irys = ^rS P , (23) 

Sp = (Q(SQJ -*Q 2 ) + {3A 2 ^5Q + A 2 {4>{&4>) - *0 2 ) + (%8Q + 4f3A 4 U5Q + A^S^j , (24) 

5q = -{QSQ + A 2 4>5(f>) : (25) 

Sp = {Q{8Q)-^Q 2 )+f3A 2 ^ 2 5Q + A 2 {^)~^^ 2 )-(^5Q + Af3A i U5Q + A A ^8^\. (26) 
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The perturbed field equations are 



'1,2 B 2 U\ { BTJ \ BU 

(5<t>y + (3H + 2(3Q)(5<t>)-+ { — + A 2 — ) 5$ = -2p [ A 2 — 6Q + (6QJ) +4*0-2^ 2 — , (27) 



dtp 
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= - 2 * ( IK + ^UP] + 4*g + 2A 2 f3^) - AA^p^Sct). (28) 
\oQ J d(p 

Integrating these equations numerically shows the perturbations of the Q-HelA are indeed suppressed. A typical 
plot of the evolution of S(f> and 8Q is shown in Fig. [T] As already said above, the difference in the behaviour of the 
two quantities is due to the large effective mass of the Q-field. Thus, the inflaton perturbations dominate and we can 
ignore the perturbations of the quintessence field. We have checked numerically that in the regime we are interested 
in, the perturbations in Q are always suppressed relative to the perturbations in <j). There is an intermediate regime, in 
which the quintessence mass is smaller but of order H and contributes to the cosmological expansion by a reasonable 
amount. In this case, the quintessence field will not sit in the minimum of the effective potential, but is attracted to 
it and its perturbations cannot be ignored. However, in this paper we do not deal with this case. 

With this in mind, the calculations proceed in the standard way, taking into account the modifications of the 
background evolution, as discussed above. The power spectrum of curvature perturbations is 



2tt J 2Air 2 M* l e \ A, 

The spectral index is found by calculating the derivative of this quantity with respect to ln(fc), where k is the 
wavenumber. Using the slow-roll condition one finds 
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FIG. 1: The evolution of the field perturbations as a function of N = In a. The upper panel shows the perturbations of the 
inflaton and the lower panel shows the perturbations of the Q-field. Here (3 = 0.5, A = 10, m = 7 x 10" 7 and k = 0.001. As 
one can see, the perturbations in Q drop to a value several orders of magnitude below that of the perturbations in <f>, due to 
the heaviness of the Q-field relative to (f>. 



The derivatives of the slow-roll parameters are: 
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So we find the spectral index to be 
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The running of the spectral index is found to be 
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For the tensor perturbations, we have 
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Writing 7-grav °^ k ng , we can show that 

n g = -2e, (37) 
whereas the tensor to scalar ratio is, using the slow-roll condition (I15p . 



7\rav 4e 2n D „ / 4/3*7 

l s I — - 



2/3/(4/3+A) 



Note that while the prediction for the scalar-tensor ratio r is modified from the standard case, the expression for the 
tensor spectral index, n g , remains the same. 



C. Consequences 

Having derived the relevant equations describing slow-roll inflation with a coupled dark energy scalar field, we will 
now discuss the consequences and predictions. The first difference from the standard case is that the expressions for 
the slow-roll parameters have changed. The origin of the modifications are two-fold: firstly, the slow-roll condition for 
the inflaton field is modified and secondly, the expansion rate is enhanced by a factor A 4 (l + 4/3/A). As a result, both 
slow-roll parameters contain an additional factor A~ 2 when compared to the standard expression, see eqns. (| 16(1 and 
(Unj) - Additionally, they are modified by a factor (f + 4/3/ A) ~ 2 (for e) and (1 + 4/3/ A) -1 (in the case of 77). Slow-roll 
inflation ends when max{e, 77} = 1. The presence of the factor of A~ 2 means that the end of inflation is delayed in this 
model, in the sense that smaller values of <fi are reached in the slow-roll phase. An immediate consequence of this is 
that the oscillations of the inflaton around its minimum - responsible for the production of particles in the reheating 
phase - will have a smaller amplitude than in the standard case. 

Another important difference is in the expressions for the cosmological perturbations. Apart from the spectral 
index of the gravitational wave power spectrum, all expressions for the perturbations have changed: the amplitude 
of scalar perturbations is different and the expressions for the spectral index n s and its running include additional 
factors which depend on the ratio /3/A. However, because the actual values of e and r\ during the last 60 efolds change 
from their values in the standard chaotic inflationary scenario, it is not obvious whether n s will be bigger or smaller 
than in the standard case. 




(a)log(m/Mpi) as a function of A and (3 (b)n s as a function of A and /3 



FIG. 2: Predictions for a chaotic inflationary potential: U = m 2 (f> 2 /2. The values of the inflaton mass (as obtained by 
normalizing the power spectrum to COBE) and spectral index vary with the parameters A and f3. In these plots, only the data 
for A > 4/3 (i.e. the region in which the inflaton is the dominant component during inflation) are plotted. 



To illustrate the effect of the dark energy field, let us consider as an example the case of a massive inflaton field 
with potential energy given by eq. |T|). The theory has three free parameters, namely m, (3 and A. We can therefore 
use eq.(|29|) and the COBE normalization to fix m as a function of (3 and A. The results are shown in Fig 



At 

the same time, we have to make sure that /3 is big enough to ensure that the Q-field sits in the minimum of the 
effective potential. We have checked numerically that for /3 > 0.05, the field sits indeed in the minimum of the effective 
potential (see fig. [3]). As one can see, the mass of the inflaton field is smaller than in the standard case (m « 10 _5 Afpi). 
This is because of the factor A e in the expression for the power spectrum and, since A > 1, a smaller mass m is 
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sufficient to obtain the correct amplitude for the power spectrum. The predictions for the spectral index n 
in Fig 2(b) As one can see, the spectral index n 



are shown 
The current 
I. Note that 



is larger than in the standard case without couplings 
constraint from WMAP alone on the spectral index is n s = 0.963^ oiti assuming a ACDM cosmology [1 
in order to compare the theory to data, one would have to study the evolution of cosmological perturbations in the 
presence of dark energy couplings in the subsequent epoch and make assumptions about the couplings to different 
matter species. The coupling (3 we have discussed so far is the coupling of dark energy to the inflaton field and is a 
priori not the same as the coupling to dark matter or neutrinos. If dark energy is subsequently only coupled to dark 
matter, then current constraints give A < 0.95 and |/3| < 0.055 at 95% confidence level Assuming that the dark 
energy coupling to the inflaton field is the same as to dark matter, this implies that the predictions of the dark energy 
coupling to the primordial perturbations are only slightly modified from the standard chaotic inflation predictions. 
On the other hand, in the case of theories with a mass-growing component [l8[, which requires larger values for A 
and P, modifications to the primordial power spectrum are predicted. In any case, when comparing the theory to 
data, assumptions about the subsequent composition of the universe and new interactions between dark energy and 
other matter forms have to be made. The results presented in this section will be useful when constructing a theory 
of inflation, dark matter and dark energy based on particle physics. 



A = 0.5 and 8 = 0.05 
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FIG. 3: These plots show the Q-field as a function N = In a (dotted lines), 
field is drawn toward its minimum (solid line). 



Even if the initial value of the field changes, the 



III. REHEATING AND PREHEATING IN THE PRESENCE OF COUPLED QUINTESSENCE 

In our scenario, when the slow-roll conditions are violated, the value of (f> is smaller than in the standard case 
and thus the oscillations of the field around its minimum are of smaller amplitude. This can affect the processes of 
reheating and preheating. In the following we will discuss in detail the effect that the presence of coupled quintessence 
has on these important epochs. 
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A. Reheating 

Any successful model of reheating relies heavily on the oscillatory behaviour of the infiaton. The oscillations are 
damped due to the expansion of the Universe by a friction term proportional to H and a term describing the decay 
of the infiaton field into radiation. Following the standard treatment of reheating (cf. (23|), the equation of motion 
for the infiaton field reads 

4>+ (3H + T T )<p+ A 2 m 2 cj)= -2/3Q<p, (39) 

where T r is the decay rate of the process. In the standard single-field case the value of the Hubble parameter is 
proportional to the square root of the total energy density, which consists of infiaton energy density and radiation 
energy density. In the present case, however, there is also the energy density of the Q-field, which contributes to the 
Hubble friction. Just at the end of inflation, for example, the ratio of the energy densities of quintessence and the 
infiaton field is roughly 4/3/A, from equation ([7]). This ratio can be potentially quite big (up to 10" 1 , or so), so the 
Q-ficld contributes significantly to the Hubble damping. This could cause heavy damping of the infiaton oscillations 
that would reduce the efficiency of the reheating process. 

In the model of elementary reheating, the equation of motion for the energy density of radiation, p r , is 

p t + AHpr = r r ( P4> + p+) , (40) 

In the standard scenario without quintessence, the energy density of the radiation produced from the decay of the 
infiaton quickly increases to a maximum value, then decreases as the decay products are diluted by the expansion 
of the universe. This continues until H f» r r , when the infiaton decays away rapidly and the radiation-dominated 
era begins with temperature, Trh- In our model, we want to avoid the dominance of the Q-field, and it is therefore 
necessary either to produce enough radiation so that [p ra d\MAX is greater than pq or to dissipate the energy of the 
Q-field quickly, so that it becomes smaller than p r . [p r ]MAX depends on the decay rate, T r and therefore on the details 
of the decay processes. For example, if the field decays to two light scalar particles with coupling gA, the highest 
decay rate is T r ~ inA [21]. This means that the amount of radiation produced is limited by the (effective) infiaton 
mass, the value of which, as was seen in Section HH must be chosen to be consistent with the COBE normalisation for 
particular values of A and (3. Figure [4] shows that if one uses constraints on the Quintessence potential from it 
is not possible to get successful reheating in this model as p r is always less than pq. 

We are thus forced to reduce the energy density of the Q-field. One method would be to let Q decay into radiation 
in a similar way to 0, with a decay rate proportional to the effective mass of the field. This would introduce a term 
p r Q, representing the energy density of the radiation produced in this manner. In this case the equations for Q and 
PrQ are: 

BV 

Q + (3H + Tq)Q + — = A 2 f3(4> 2 - AA 2 U{<j>)), (41) 
p rQ + 4H PlQ = T Q {p Q + PQ ) (42) 

However, numerical calculations (see fig. [5} show that the radiation produced is still not sufficient. This is partly due 
to the small magnitude of the term driving the production of radiation: Tq (pq + pq) — TqQ 2 . (Q is small as the 
effective potential is very flat at this stage.) 

A favourable alternative is to vary the parameter A that controls the steepness of the potential of the Q-field. A 
larger A means that the field rolls more quickly to larger values. This has a twofold effect: [p ra d\MAX is increased 
(as the damping is reduced, so the source term in eqn. (|40[) does not decrease as quickly as in the previous case) 
and pq is decreased, so [p ra d\MAX does not have to be very large for radiation to dominate the universe. It can be 
seen in fig. [6] that using a large value of A does indeed allow radiation to dominate the universe. It is interesting 
to note that this is true even if T r is much less than the maximum possible decay rate. In this case, as can be seen 
in the lower panel of fig. [(5J the total amount of radiation produced is dramatically reduced, leading to a smaller 
reheating temperature than in the case of large T r . The same result is true in the standard single field case, where 
T RH ~ O.55.g7 1/4 (V8^Mp;r0) 1/2 23]. The difference is that in our model if the decay rate is too small, one cannot 
satisfy the condition p r > pq. 

It is clear that if our model is to be consistent, the parameters of the quintessence field must be constrained by 
this theoretical consideration and A must be large. It has been shown [18| that it is possible to achieve a realistic 
cosmology with large A in the case with a growing matter component. As Figure [S] shows, it is possible for radiation 
to dominate the universe after the decay of the infiaton using values that satisfy this lower bound. 
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A = 0.5 case 
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FIG. 4: The log of the ratio of the radiation energy density to the total energy density for different values of A and /3. It is 
not possible to get radiation dominance using constraints on the parameters from [l9|. We have restricted ourselves to the case 
where the coupling, /3, is large enough make Q stay in the minimum of its effective potential. 



B. Preheating 

We have seen that reheating is considerably affected by the presence of a coupled quintessence field. In particular, 
the process is less efficient than in the standard case, because the Hubble damping is much larger. In the following, 
we discuss the process of preheating, in which the inflaton energy is converted into other particles in a very efficient 
way. We consider the possibility that the energy density of these decay products could counterbalance the effect of 
the Q-dominance and lead to successful reheating. 

Following the standard description of preheating, as discussed in [24|, [2^, [26[, we assume that the inflaton field 
decays to a light scalar field x °f mass m x . To describe the interaction between and x, we will add the following 
term in the Lagrangian: 

A£ int = - 5 A0 X 2 . (43) 

The conformal factor A arises because our analysis is performed in the Einstein frame. This term describes a three- 
legged interaction involving one inflaton particle and two x particles. Another common choice of interaction is of the 
form —^g 2 (p 2 x 2 - As similar behaviour arises in both cases and we are interested in the more general consequences 
of the presence of the Q-ficld on the mechanism of reheating, we will concentrate our attention on the interaction in 
eqn. (|4"3"]) . We assume that the effective coupling constant gA is less than the frequency of the inflaton oscillations 
mA, so that the interaction is not modified too much by quantum corrections. Thus, we have the limit g < m. The 
vacuum expectation value of the x-field is zero, so the Friedmann equation and the classical equation of motion for 
the inflaton will be unaffected. Expanding x around zero, we find the following equation for the perturbations of x, 
which are interpreted as particles after quantization: 
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FIG. 5: In this simulation, the model was modified to allow the Q-field to decay to radiation in a manner analogous to 4> by 
adding a friction term Tq ~ m c g to eqn. (|3]). The radiation produced in this manner contributes a fraction of the total energy 
density. These plots shows the log of the fraction of the energy density accounted for by the total amount of radiation (i.e. 
including also that from the decay of the inflaton). 
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Xk + 'iHxk + (K + ™\ + ^gA^j X k = 0. (44) 

To study the resonance, we neglect the expansion of space and introduce a sinusoidal ansatz with which to describe 
the oscillations of the inflaton, <f> = &sm(mAt). In this case, eqn. (|44p can be written in the form of a Mathieu 
equation, i.e. 

x£+[A fc -2«cos(2z)]xfc=0, (45) 

with 



4(fc 2 + 
m 2 A 2 



Ak = (46) 



1 = — 27' ( 4? ) 
m A 

z = mAt/2, (48) 

where prime indicates a derivative with respect to z. These are the standard equations with the addition of the 
factors A in the expressions for Ak, q and z (cf. [2l[ or [ll])- As we have also seen in Section HH the inflaton mass is 
smaller than in the standard theory once the theory is normalized to COBE and the amplitude $ of the inflaton field 
is smaller. We note that for a given value of g, the parameter q, which determines the behaviour of the solutions to 
the equation above, can be considerably bigger than in the standard case. 
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FIG. 6: These plots show the evolution of p^, pQ and p r (on a log scale) for large A. Here, A = 10, /3 = 0.5 and m = 7 x 10 -7 . 
In the upper panel, the decay rate is large (T r — 0.9mA) so radiation domination occurs very quickly. In the lower panel, 
although the decay rate is much smaller (T r = 0.001mA), the radiation can still dominate over the Q-field but this process is 
of longer duration. The reheating temperature would be smaller in this latter case. 



A similar equation can be derived for the Q-particles produced, since the Q-field couples to the inflaton field as 
well. We obtain (Qk is the Fourier component of the perturbation of the Q-field around its VEV) 

Ql + [4 Q) - 2g«> cos(2* Q )]Q & = 0, (49) 

with 

4 Q) = ^# + >, (50) 
k m 2 A 2 5 

q iQ) = ^ 2 f3 2 A 2 , (51) 
z Q = mAt. (52) 

Since tuq is much bigger than mi, we see that A]l is dominated by the first term and therefore A^ ^> q(®\ 
This means that the periodic term in the equation for Qk does not play an important role and there is no resonance 
production of Q particles. Therefore we concentrate on the production of the light \ particles in the following. 

Let us consider first the case of narrow resonance, in which g« 1. Solutions to eq. (f45|) falling within particular 
resonance bands are exponentially unstable and take the form \k ~ exp(/ifcz). The most important band is the first 
one, for which the resonance reaches its maximum at fi^ ~ The number density of particles with momentum k is 

f\x\ 2 , | , 2 \ 1 ,„v 

nfe = T Ur + lxl rr (53) 



where ui 2 , = A^ — 2qcos(2z). From this, one can see that, while narrow resonance continues, the number of particles 
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with k ss mA/2 increases exponentially, with n& ~ exp(gz) (as x is, by design, much lighter than the inflaton and the 
resonance occurs at ~ 1). 

Narrow resonance will only be an important decay mechanism as long as it is more efficient than the perturbative 
decay, which corresponds to the elementary theory of reheating discussed in the last section. In narrow resonance, 
the number of particles in the resonance band k increases exponentially as rih ~ exp(gz) ~ exp(qmAt/2) so the decay 
rate can be approximated by qmA. In the regime where perturbative decay itself is efficient (i.e. T x > H), we find 
that the condition for narrow resonance to be the leading effect is 

qmA > T x , (54) 

where T x = g 2 A/ (Sum) is the decay rate as obtained by standard quantum field theory methods. The second 
phenomenon that affects the timescale of narrow resonance is the redshift of momenta out of the resonance bands. 
The exponential nature of the resonance means that the rate of production of x particles depends on the number of 
particles present already. If the modes do not spend enough time within the resonance band, rik will be small and 
the process will be inefficient. 

The width of the resonance band at k = mA/2a is kq. If we assume that the resonance is most efficient in the 
middle of this band we have Ak = mAq/2 . We wish to find the time, At, that the mode spends in the region 
Afc. Using Afc = |4||At = kHAt, this is At = qH~ x . During this interval, the number of particles increases as 
e-xjp(qmAAt/2) = exp(q 2 mA/2H). So the second condition required for decay by narrow resonance to be efficient is 

q 2 mA > H. (55) 

We can write the inequalities in terms of the inflaton amplitude, <I>. 

<*> > (56) 



m V mAH 

$ > • 57 

4.9 

A small value of m means that parametric resonance can continue unhindered for long enough to produce large 
amounts of x particles, which subsequently decay to radiation. 

If, instead, the parameter q is greater than one, we are in the regime of broad resonance. Here, the mass term in 
eqn. (|4"5|) is dominated by the sinusoidal behaviour of the inflaton. When this term is zero, the number of particles 
given by eq. ([53]) is not conserved and one observes resonance in a broad range of modes. This can be seen in figure 
[7J The small value of m required to match the COBE normalisation means that the parameter q is very large in this 
case, leading to a large increase in the number of x particles, despite the small amplitude of the oscillations. 

The theory governing the energy density of the radiation produced by this method and the backreaction of the 
particles on the evolution of the inflaton is dependent on the model with which one is working and the free parameters 
therein. We can assume, at the very least, that the energy density of the products of parametric resonance will not 

exceed that of the inflaton prior to decay. As, at the start of the reheating era, pq « P4>, Hubble damping of the 

inflaton can reduce the energy density of the inflaton field to a value less than that of the Q-field before the field can 
complete many oscillations about zero. This rapid damping is exhibited in figure [71 Because any model of parametric 
resonance assumes the inflaton is oscillating, by the time a significant quantity of x-p&rticles are produced, the inflaton 
is incapable of producing enough radiation for it to dominate over the Q-field. Therefore, even if we employ parametric 
resonance, we still have to use a large value of A to reheat the universe in our model of coupled quintessence. 



IV. CONCLUSIONS 



We have described a model of coupled quintessence in which the scalar field responsible for dark energy is coupled 
to the inflaton field. To be specific, we have used the exponential potential in eqn. and assumed that the coupling 
between the fields is of the form A(Q) — exp(f3Q/Mpi) with f3 positive, so that the effective potential has a minimum. 
The large effective mass of the Q-field drives it into this minimum for a large range of A and (3 values, and in Section 
III Al we used this to evaluate the slow-roll parameters. We found that the expressions for rj and e are modified by 
factors of A~ 2 (l + 4/3/A)~ 1 and A~ 2 (l + 4(3/\)~ 2 respectively, meaning that the end of inflation is delayed in this 
model, so that the inflaton field has a smaller field amplitude during its oscillations. 

We further discussed the modification to the cosmological perturbations. The equations for the spectral index and 
its running (in terms of the new slow-roll parameters) vary with the ratio /3/A. Normalising the power spectrum to 
COBE, we found that the mass of the inflaton field is required to be smaller than in the standard case, whilst the 
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FIG. 7: The oscillations of the <j> field are still capable of driving parametric resonance, despite their small amplitude. In 
this plot, f3 = 0.1 and A = 1. g is set to 0.9m, where m is the inflaton mass, which is chosen to be 4 x 10 -8 to match the 
COBE normalisation. In upper plot we show the behaviour of <j>, whereas the lower plot shows the number density rik of the x 
particles. The dotted lines on the upper plot show the bounds given in eqns. (|56p and (|57p . below which narrow resonance is 
not possible. 

effect of the presence of a dark energy coupling is to increase the value of n s , for a given A. The parameters f3 and 
A can be constrained by observations if one assumes that the coupling between the inflaton and dark energy has the 
same equal magnitude as the coupling between dark matter (or neutrinos) and dark energy. 

We found that the presence of the Q field during reheating leads to an unnaturally large amount of Hubble damping 
in the equation for the inflaton oscillations during reheating. We showed that if reheating is to be successful, A must be 
large. The dark energy field rolls faster down its potential, which in turn reduces the Hubble damping and allows the 
radiation produced during reheating to become the dominant component of the universe. We also considered how the 
mechanism of parametric resonance during preheating might work in this model. We found that parametric resonance 
could still produce large amounts of radiation as the smaller value of m required to match the COBE normalisation 
decreases the lower bound on the amplitude of the inflaton oscillations required for resonance. However, if A is small, 
the energy density of the inflaton quickly becomes less than that of the Q-field so a large value of A is still required 
to reheat the universe in our model. The fact that a large A is required for successful reheating has interesting 
consequences. For example, it is not compatible with the condition A < 0.95 stated in (l9| . Such models would be 
only consistent if the quintessence field does not couple to the inflaton field. Alternatively, the dark energy field could 
couple to a subdominant species, such as neutrinos, as e.g. in [l8l |. 
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